Fingerprints of Primordial Universe Paradigms 
as Features in Density Perturbations 



Xingang Chen 

Center for Theoretical Cosmology, 
Department of Applied Mathematics and Theoretical Physics, 
University of Cambridge, Cambridge CBS OWA, UK 

Abstract 

Experimentally distinguishing different primordial universe paradigms that lead to the Big 
Bang model is an outstanding challenge in modern cosmology and astrophysics. We show 
that a generic type of signals that exist in primordial universe models can be used for such 
purpose. These signals are induced by tiny oscillations of massive fields and manifest as fea- 
tures in primordial density perturbations. They are capable of recording the time-dependence 
of the scale factor of the primordial universe, and therefore provide direct evidence for specific 
paradigm. 



Current observations have mapped out a detailed 13.7 billion years of Big Bang expansion 
history for our universe. One important knowledge we learned from these observations is 
that, at the very beginning of the Big Bang, there are tiny fluctuations in an otherwise very 
homogeneous and isotropic background density. These fluctuations, called the primordial 
density perturbations, become the seeds for the subsequent evolution of the large scale 
structures. They turn out to have very special properties. They are seeded at super-horizon 
scales, and are approximately scale- invariant, Gaussian and adiabatic [T]. Understanding the 
origin of these initial conditions and hence the origin of the Big Bang becomes an outstanding 
challenge for the modern cosmology. 

The inflation paradigm prH] is the leading candidate for creating such a primordial uni- 
verse. The simplest inflationary scenarios not only explain why our universe is homogeneous 
and isotropic, but also the properties of the density perturbations. Nonetheless, based on 
current data, there exist important ambiguities and degeneracies. Competing paradigms 
have been proposed to produce the same initial conditions. These alternatives include the 
paradigms of cyclic universe [5J, matter bounce [HIE] and string gas [8]. Arguably, the infla- 
tion still remains as the best available paradigm, because its generic predictions naturally fit 
the data and its microscopic origin in terms of fundamental physics is promising. However, 
such opinions may be subject to personal tastes; they are mo del- dependent and may even 
evolve with time. Besides theoretical establishment for each paradigm, a universal and un- 
ambiguous distinguisher is in terms of data. Namely we need to find some properties in the 
density perturbations that can clearly distinguish different paradigms. Such distinguishers 
are desirable even if one believes in a paradigm. In order to do so, such properties should 
satisfy the following requirements. They have to be shared by all general models in one 
paradigm, not just by a small subset of models; and they have to be distinctive for different 
paradigms. 

So far the primordial tensor mode is regarded as the only possible solution to this ques- 
tion. They originate from gravitational waves in the primordial universe [9], and may be 
observed in terms of polarizations in the cosmic microwave background (CMB) [TUtfTT]. 
Inflation models indeed have general predictions on the tensor mode. Namely, the tensor 
modes are approximately scale-invariant with a red tilt; for some models they are observable. 
Some alternative paradigms predict non-observable tensor modes, such as the cyclic universe 
paradigm [5J, or observable ones with blue-tilt, such as the string gas paradigm [12] . However 
there exist important caveats. Firstly, even for inflation, tensor modes sourced by the pri- 
mordial gravitational waves are not guaranteed to be observable. While the best sensitivity 
for the tensor-to-scalar ratio achievable by experiments in the near future is Ar ~ 0(l(r 3 ), 
the inflation models predict anywhere between r ~ (9(10 -1 ), for large field models, and 
r ~ (9(10~ 55 ), for small field models with TeV-scale reheating energy. Secondly, if we con- 
sider more general alternatives, scale-invariant and observable tensor modes are possible. 
The equation of motion obeyed by each polarization component of the tensor modes is the 
same as that by the massless scalar with unit sound speed. So the tensor modes can be scale- 
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invariant even in non-inflationary spacetime, just as the scalar. They become observable if 
the Hubble parameter is large. The matter contraction paradigm is such an example [Sl[7]. 

Therefore it is very important to look for complimentary properties in the density pertur- 
bations that can be used as a model-independent distinguisher between different paradigms. 
This is the main purpose of this paper. 

Almost by definition, different primordial universe paradigms are distinctively character- 
ized by different time-dependence in the scale factor a(t). To be explicit, we approximate 
the arbitrary time-dependent scale factor as the general power-law, 



The density perturbations are seeded at superhorizon scales, so we require that the quantum 
fluctuations exit the event horizon \ar\, where r is the conformal time defined by dt = adr 
and is related to t by ar = t/(l — p). To satisfy this requirement, for p > 1 we need an 
expansion phase, so t runs from to +00; for < p < 1 we need a contraction phase, so t 
runs from —00 to 0; for p < 0, we again need an expansion phase, so t runs from —00 to 0. 
For example, p > 1 corresponds to the inflation, p = 2/3 the matter contraction, < p <C 1 
the ekpyrotic (slowly contracting) phase, and -1 C p < the slowly expanding phase, r 
always runs from —00 to 0. 

Many properties of the density perturbations are convoluted consequence of this scale 
factor, and this is a primary reason for the ambiguity. Perhaps the most direct way to 
distinguish different paradigms would be to find some properties that can directly record the 
time-dependence of the scale factor. In this paper, we show that a type of feature models 
that generically exist in the primordial universe paradigms provide such opportunities. 

We first describe this idea heuristically. Consider models with some small and repeated 
features. Such features can be various kinds of small structures in the realistic primordial 
universe models, and we will give a very general example shortly. These features generate 
small and repeated beats in cosmological parameters during the evolution, which can be used 
as a clock and related to the time t. In the mean while, the scale factor a and the comoving 
momenta k of quantum fluctuations are related by the physical momenta by definition. The 
way that the features imprint themselves on the density perturbations is through the beats 
on the physical modes. Therefore the scale factor a as a function of t is translated to the 
comoving momenta k as a function of features. What we can measure is features as a function 
of k. Knowing that we know the inverse function of a(t). Let us now examine how this idea 
may be implemented and how general it is. 

In order to achieve the goal, we need to look for standard physical clocks. Such clocks 
should generate repeated features with known time-dependence, although not necessarily 
periodic. They should exist generically so the signals are more likely to be observed, and 
they should also be associated with a set of specific patterns identifiable in observations. 
The classical oscillation of massive fields are ideal for such purposes. 
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Massive fields with mass much larger than the horizon mass-scale l/|ar| exist ubiquitously 
in models of primordial universe. Even when we think of single field models, in the context of 
a unification theory, what we really have in mind is models with many massive modes. The 
single field model is obtained as the low energy limit at or below the energy scale l/|ar|, after 
the massive modes are integrated out. If we take into account these massive field directions, 
the trajectory of the effective single field in this multi-dimensional field space is unlikely 
a straight line. Generically it should turn from time to time. During the turning, the field 
4> and the massive modes couple, and this induces the energy transfer. The massive fields 
are excited and oscillate around the potential minima classically. Because generically the 
fraction of energy transfer is small, the oscillation amplitudes are small but the frequency 
are large. So for most purposes, these effect are averaged out and safely ignored. But for 
our purpose, these small oscillations are a good candidate for the physical clock that we 
are looking for. So let us examine more closely the behavior of the massive field a in the 
power-law background. 

The equation of motion is 



where the Hubble parameter H = p/t. The solution is given in terms of Bessel functions. 
The asymptotic behavior of these Bessel functions at the limit m a t ^> p 2 is given in terms 
of sinusoidal functions, and we use these to approximate the oscillatory behavior of a, 



where a is a phase, and a a is the initial oscillation amplitude at t — to- Such oscillations 
induce an oscillatory component to the Hubble parameter H, because 



The leading term on the right hand side of (JH) does not oscillate in time because the energy 
is converting between kinetic and potential energy back and forth and conserved in the 
leading order. The oscillatory component for H, which we denote as H osci , comes from the 
subleading terms. Using ([3]), we get 



a + 3H& + m\o = , 



(2) 




(3) 
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This in turn induces the oscillatory components for the parameters e = —H/H 2 and rj 
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e/(He). Again we use the subscript "osci" to denote their oscillatory components, 

2 2 / j. \ — 3r> 

a\mi ( t x 



- cos(2m a t + 2a) , (6) 



^osd = -]^^3^J cos(2m CT t + 2a) . (7) 

Generally, multiple massive fields are excited sporadically along the trajectory with dif- 
ferent spectra. These oscillations induce corresponding beats in the cosmological parameters. 
The question now is how sensitive the observables are to these oscillations and how they can 
be identified in experiments. To study this, let us turn to the subject of Bunch-Davies (BD) 
vacuum and resonance mechanism. 

In nearly all models of primordial universe, the quantum fluctuations start their life in a 
vacuum that is mostly BD. These fluctuations later exit the event horizon and become the 
seeds for the large scale structure. Their general property can be seen as follows. Consider 
the fluctuations of an effectively massless scalar field, 50(x, t), in a general time- dependent 
background with the scale factor a(t), 

£ = y(^) 2 - \{dM? ■ (8) 

A quantum fluctuation with comoving momentum k is within the event horizon if k > l/|r|. 
In this limit, the equation of motion for the fluctuations approaches that in the Minkowski 
spacetime limit. Along with the quantization condition, 

a 3 5<f)5<j)* - c.c. = i , (9) 

the mode function in the subhorizon limit becomes 

5(f) —^=e~ ikT . (10) 
aV2k 



We have chosen the positive-energy mode, which corresponds to the ground state of the 
Minkowski spacetime, to be the BD vacuum. The effect of the background time- dependence 
is incorporated adiabatically in fflQ|) . The most important and universal property of fflQ|) is 
the oscillatory factor e~ lkT . 

While the superhorizon evolution of different paradigms are very mo del- dependent, the 
subhorizon BD-vacuum behavior is universal. This applies to both inflationary and non- 
inflationary scenario, expansion and contraction universe, attractor and non-attractor evo- 
lution, single field and multifield model, and curvaton and isocurvaton modes. 

If the background cosmological parameters have small oscillatory components with fre- 
quencies much larger than the energy scale of the event-horizon, they introduce a high 
energy probe to the BD vacuum by resonating with the vacuum component that has the 
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same frequency [IB] . Since the BD vacuum has the background time-dependence, different 
momentum modes get resonated at different time. This process produces large correlation 
amplitudes and distinctive scale-dependence (running). They manifest as features with spe- 
cific patterns in density perturbations, which we will compute shortly. 

To summarize, we see that three universal properties nicely fit into each other for our 
purpose. Firstly, the generic oscillation of fields with mass much larger than 1/|ot| can 
be used as the standard clock, and they generate small beats with known time-dependence 
in cosmological parameters. Secondly, these oscillations affect the density perturbations 
through the universal BD-vacuum instead of the model-dependent superhorizon behavior, 
which makes general analyses possible. Thirdly, the sub horizon scale is precisely where the 
strong resonance effect takes place, generating large effect on density perturbations. 

We consider a two-field model in the general time- dependent background as the explicit 
example. The field is responsible for the leading order scale-invariant power spectrum 
P^o- This field may be the one driving the background evolution, or a different field that is 
introduced to source the density perturbations only. The details are highly model-dependent, 
but this is not our concern here. We are interested in computing the resonance effects induced 
by massive fields as corrections to the leading power spectrum and as the leading bispectra. 
To be as general as possible, after the massive field a is excited and starts to oscillate at 
time t , we consider the massive field as a spectator which couples to only through gravity, 



More complicated couplings are of course possible, but this is the minimum case. The 
correlation functions for the scalar perturbation £ are computed in the in-in formalism. See 
Ref. [H] for computational details and Ref. [15] for general methods. For power spectrum, 
we have 
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For bispectrum, we have 
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We will quote the bispectrum in terms of a function S(ki, k 2 , k 3 ) defined as [15] 



(CkiCkaCO = S(k u k 2 ,k 3 ) 



1 



-P c 2 (2tt) 7 5 3 £>) . 
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The mode function u^(t) in these expressions is the Fourier transform of £(x, r), and as 
mentioned only the BD behavior is needed in our computation. In the above integrals, the 
resonance happens between the oscillatory components of the parameters, e osci and rj OSC i, and 
the mode functions. 

For arbitrary p, both power spectra and bispectra have similar p-dependent resonant 
running behavior, 
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The resonant running refers to the scale-dependence of the functions sin[. . . ] in these ex- 
pressions. We have also included the typical scale-dependence of the amplitudes in front of 
these functions. The parameter k r denotes the first resonating .fT-mode when the massive 
field starts to oscillate (K = 2ki for power spectrum and K = k\ + ki + k 3 for bispectrum). 
For the expanding background, p > 1 and p < 0, lower X-modes resonate earlier and the 
above formulae apply to K > k r ; for the contracting background, < p < 1, larger fc-modes 
resonate earlier and they apply to K < k r . 

The most important property in the above results is that the resonant running for different 
paradigms are distinctive. The differences are clear even within a couple of efolds. This 
demonstrates the idea of using features to probe the scale factor. The arguments of the 
sinusoidal functions are power-law functions of comoving momenta with a power 1/p, which 
is the inverse function of the power-law in the scale factor. 

Generally speaking the overall sizes of these signals depend on detailed models. If we 
restrict to the inflation case, p>l, the full results become 
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(18) 



where the phase a = (2m a / H)(l — ln2m CT ) + 2a + it /A. We will use these to demonstrate 
the sizes of the signals in the inflation case. The resonant running in this limit reproduces 
the form found by Chen, Easther and Lim (CEL) for inflation with periodic features p2]. In 
retrospect, the fact that the arguments are proportional to ln.fr tells us that this is inflation. 

At least for inflation models, the amplitudes of the resonant power spectra and bispectra, 
which we label as (AP(/P( )a and f^L respectively, can be easily made very large. We denote 
the fraction of the kinetic energy of <j), that is converted to the energy in the cx-field during 
the turning and induces its oscillation, as 0. So m^a^ ~ (3<fi 2 . Also note e ~ <fi 2 /(MpH 2 ). 
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From (I17p and fll8p . we have 



Even for a tiny fraction of energy transfer, the resonance amplitudes can be quite large. For 
example, for (3 ~ 1CT 2 , m u /H ~ 10 2 , we have AP^/P^ ~ 0.1 and f^L ~ 10 3 . The resonance 
mechanism plays a crucial role in producing the large amplitudes. With same (3, increasing 
the mass m CT reduces the oscillation amplitude a a, but this is overcome by the strength of 
the resonance which is proportional to higher powers of m a . 

It is important to note that the distinctive oscillatory running behavior in the above 
resonant forms will not be changed by curvaton-isocurvaton couplings in multifield evolution. 
Any effect that oscillates faster than the horizon time-scale generates additional resonance 
forms that superimpose onto each other. Any effect that varies much slower can only change 
the overall envelop of the resonance forms, by either changing their overall sizes or introducing 
scale-dependent modulations. This is another crucial point that makes this type of features 
a model-independent general distinguisher between paradigms. 

There are also interesting scale-dependence in the amplitudes. These dependence is much 
milder in the sense that the range of scales AK over which significant variation takes places 
is larger than or comparable to the local scale K itself. For general multifield models, such 
scale dependence can be changed due to scale-dependence in curvaton-isocurvaton couplings. 
So they may not be as faithful as the resonant running to the specific paradigm. Nonetheless, 
as we will see, there are some very sharp differences between paradigms that can be used as 
supportive evidences. 

It is also important to note the caveats in our arguments. In principle one may be able 
to engineer artificial features in a model that mimic the signals we get so far. The type of 
clocks we considered are periodic; but one may engineer non-periodic background oscillation 
that resonates in a non-inflationary background, and they conspire to give the same resonant 
running of the CEL form as in (1181) . For instance, the oscillation of the form ~ e *9 ln (*/*o) j n 
the non- inflationary power-law background induces the same CEL form. Another possible 
caveat is that the mass may be time-dependent so the period of the clock changes. 

These ambiguities can be reduced by looking for multiple characteristic signals. The main 
point is to strength the evidence that these features are induced by periodically oscillating 
massive fields. We have so far only concentrated on the resonant running behavior. There 
are at least two types of multiple signals characteristic of massive fields. 

Firstly, the density of massive fields is proportional to a -3 ; in addition, they induce a con- 
stant resonance scale while the event-horizon sizes have different time-dependence for differ- 
ent paradigms. These properties determine the different running behavior of the amplitudes 
in (1151) and (1161) . In particular, inflationary background has a very red running index —3, 
while the Ekpyrosis has a very blue running index oc 1/p. Although these types of running 
are less robust than the resonant running, this dramatic difference can be used as supportive 
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Figure 1: The first type of fingerprints for different paradigms: multiple oscillations. From 
top to bottom: p = 10 (inflation), p = 2/3 (matter contraction), p = 0.2 (Ekpyrosis), 
p = —0.2 (slow expansion). The red/dark spectra correspond to m a /H = 10, 20, 120 excited 
by a common sharp feature; the green/ligh spectra correspond to m a /H = 30, 150 excited 
by another common sharp feature. 



evidence. These running amplitudes will be modified if the mass is time-dependent. For 
example, for inflation, to make the resonant running deviate significantly from the CEL 
form, we need m a / (m a H) > 0(1); but this introduces dramatic scale-dependence in ( |T9l . 
In general, since the massive fields are excited sporadically at different places with different 
spectra, we expect that each paradigm has its unique fingerprints consisting of such multiple 
signals of the same type. This is illustrated in Fig. [TJ 

Secondly, oscillations of the massive fields are triggered by some sharp physical pro- 
cesses, such as a sharp turning (with large angular velocity 9 ^> H but not necessary large 
turning angle). Such sharp features are also associated with characteristic observational 
signatures. In contrast to the resonance case, sharp feature signals are universal to different 
time- dependent backgrounds. See Ref. [14] for details. The characteristic property is the 
sinusoidal running 

~ sm(K/k + phase) , (20) 

where ko is the starting point of this running and it is related to the conformal time tq 
at the location of the feature by k = |r | _1 . Such signals cannot be used to distinguish 



8 




Figure 2: The second type of fingerprints for different paradigms: in conjunction with sharp 
feature. From top to bottom: p = 10 (inflation), p = 2/3 (matter contraction), p = 0.2 
(Ekpyrosis), p = —0.2 (slow expansion). The Green/light spectra are generated by a sharp 
feature at K = ko = 100 and have sinusoidal running; the blue/dark spectra correspond to 
two massive fields (m CT = 25, 60) excited by this sharp feature and have resonant running. 

different paradigms, but they can be used to identify the location of the sharp feature and 
the subsequent massive modes oscillations. Because of this connection, the parameters for 
these two different runnings, (TIE]) and ( 1201) . are related, 

k r _ p 2m a 

k ~\l-p\ H ' [ZL) 

where H Q is the Hubble parameter evaluated at r . This is illustrated in Fig. [2j 

These signals present special opportunities and challenges for experiments and data anal- 
yses. Typically they are highly oscillatory in i^-space, which means that in the multipole 
space the binning interval A£ has to be much smaller than the local £, A£/£ ~ H/m a <C 1. So 
experiments capable of measuring density perturbations in high multipoles are required. The 
Planck satellite is observing the CMB at multipoles of a few thousands. The ground-based 
telescopes, South Pole Telescope (SPT) p3] and Atacama Cosmology Telescope (ACT) [17] , 
can go up to ten thousands. More speculatively, the 21cm hydrogen lines may be observed 
in much lower redshift and with much higher multipoles. 

The CEL form, and most other resonance forms, have highly oscillatory and character- 
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istic running behavior. This can easily make them orthogonal to many other primordial 
signals and late-time contaminants. Non-linear gravity effects and CMB evolution limit the 
detection sensitivity for conventional primordial bispectra to be Jnl ~ 1. For SPT and 
ACT, astrophysical effects such as point sources and Sunyaev-Zeldovich effect dominate over 
the primordial information for i > 2500. However, since the resonant forms are orthogonal 
to these non-primordial contaminants, new assessment is necessary to find out the main lim- 
iting factors and make forecasts; and new observationally accessible windows can be opened 
up to probe the primordial universe even for existing experiments. 

Conventional data analyses methods start with a specific separable template. But here a 
scan over families of non-separable functional forms are necessary. The modal decomposition 
method developed by Fergusson, Shellard and Liguori [T8lfT9] is well-suited for such purposes. 
So far it has been mainly applied to cases with less dramatic scale dependence, but should 
be able to be generalized to cases relevant to the signals studied here. 
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